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Abstract 


We propose a numerical scheme to solve a general class of time-fractional 
order telegraph equation in multidimensions using collocation points nodes 
and approximating the solution using double shifted Jacobi polynomials. The 
main characteristic behind this approach is to investigate a time-space collo- 
cation approximation for temporal and spatial discretizations. The applica- 
bility and accuracy of the present technique have been examined by the given 
numerical examples in this paper. By means of these numerical examples, 
we ensure that the present technique is simple, applicable, and accurate. 
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1 Introduction 


Fractional differential equations [2, 20] are exhibited as capable mathemati- 
cal tools for factual and more precise depiction of various phenomena. They 
show up in different territories, counting mathematical chemistry [12, 19], 
viscoelasticity [25], biology [20], electrochemistry, physics [17], semiconduc- 
tors, seismology, scattering theory, heat conduction, fluid flow, metallurgy, 
population dynamics, optimal control theory, mathematical economics, and 
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chemical reaction. As the increasing of employing fractional partial differen- 
tial equations [21, 10, 30] in many social and scientific fields, the principle 
challenge we defy is that getting answers for them. Unfortunately, for most 
of these fractional partial differential equations, no one able to achieve an- 
alytic solutions for such problems. There are an extraordinary number of 
demonstrating and fractional order differential equations, which have been 
illuminated numerically utilizing different methods, see [5, 6, 1, 8, 23, 24]. 

There are many applications of the telegraph equation such as signal anal- 
ysis for transmission, propagation of electrical signals, and modeling of the 
reaction diffusion. Sinc-Legendre collocation method [26] has been used to 
solve the time-fractional order telegraph equation (T-FOTE). By means of 
radial basis functions, Hosseini, Chen, and Avazzadeh [18] treated with time- 
fractional telegraph equation. Wei et al. [27] applied the fully discrete local 
discontinuous Galerkin method to solve fractional telegraph equation. In [28], 
space- and time-fractional telegraph equations have been solved by using ho- 
motopy perturbation method. Furthermore, the authors in [16] transformed 
a wavelet method based on Haar wavelets to solve space and time fractional 
telegraph equations. 

In this work, we mean to build up some successful and productive collo- 
cation schemes to comprehend time-fractional order telegraph equation. One 
great advantage of such schemes is that it reduces the problems under con- 
federation to systems of algebraic equations by using combination of basis 
functions of shifted Jacobi polynomials and the Gauss-shifted Jacobi nodes 
as the collocation nodes. The collocation method has successfully been ap- 
plied to many situations [9, 4, 3, 14, 13, 15]. The main advantage of the 
proposed method is that is easy to implement, and also, we obtain highly 
accurate semi-analytic solutions via few number of retained modes. 

The framework of this paper is as per the following: In the following area, 
we present few relevant properties of fractional derivatives and shifted Jacobi 
polynomials in the coming section. Section 3 is assigned to the theoretical 
derivation of the shifted Jacobi collocation (SJC) method for one-dimensional 
T-FOTE with the homogencous and nonhomogeneous conditions. Section 4 
is assigned to applying the SJC method for two-dimensional T-FOTE with 
the homogeneous and nonhomogeneous conditions. Moreover, in Section 6, 
several numerical examples and simulations are presented to clarify the effec- 
tiveness and accuracy of the proposed underlying method. Finally, related 
conclusions and observations are introduced. 


2 Mathematical preliminaries 


Some definitions and preliminaries related to fractional calculus [25] are 
stated in this section. Also, we listed some properties related to the shifted 
Jacobi polynomials. 
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2.1 The fractional integration 


For fractional integration of order ~ > 0, we can find different definitions, 
which are not necessarily equivalent; see [22]. The more used definitions are 
Riemann-—Liouville and Caputo fractional definitions. 


Definition 1. The Riemann—Liouville integral of order > 0 is defined 
as 


J" f(a) = Ww [ Gash GS) est, 


J° f(x) = f(z), 


The operator J" satisfies the following properties: 


(1) 


JET’ f(x) =IU f(a), INI’ f(a) = IYI" f(a), 
itl) pep, (2) 


Jtg? = —s_—_ ; 
“ TB+14p) 


Definition 2. The Riemann—Liouville fractional derivatives of order vy is 
obtained by 


1 qd™ 
T(m — p) dz™ 


D¥ f(x) = fw ty™-/-1 Fe)dt), m-—l<p<m, «>0, 
0 


here m is used as the ceiling function of w. 


Definition 3. The Caputo fractional derivatives of order yu is defined as 


1 as a™ 
Cpe m—p-l = = 
D¥ f(x) = <6 oo (a —t) ere m—-1l<u<m, c>0. 
With simple calculations, we obtain 
D'c =0, (C is a constant) (3) 
0, form € No and m < [yu], 
D'e™ = T'(m + 1) hc form € No and m > [pu] (4) 
T(m+1-— yp) orm ¢ N andm > |], 


where [4] and || are used as usual for the ceiling and floor functions, 
respectively. 
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2.2 Properties of shifted Jacobi polynomials 


By means of the main properties of Jacobi polynomials, we conclude the 


following: 


POG) =i. PPP a= sas 64 ata p), 


(-1)*T(k +6 +1) (5) 


PRP a) = (Pea), PP(-1) = “a 


where a, 6 > —1, a € [-1,1], and 


(08) o (Qkt+at+8+1)(2k+a+ 642) 


oe 2(k+1)(kK+a+6+1) 

p(or8) (6? - 0?)(2k+a+ 8 +1) 

RO  8(k+D(kK+a+6+1)(2k+ a4 8)’ 
(os) _ (kt al(k+ 6)(2k+at+B+2) 
. (kK+1)(kK+a+6+1)(2k+ a+ B) 


Furthermore, the rth derivative of a (x), is computed as 


Tjt+at+6t+qt 1) B(a-+r.8-+r) (x) (6) 


D’P(o#) = : 
= otG+a+B+l) 2 


J 


where r is an integer. For the shifted Jacobi polynomial pie (x) = 


poe —1), £> 0, the explicit analytic form is written as 


k : 
(a8) _ pag TEBE OE ROE +1) 
Pek @)= LCV TG arirerats+De-ye™ 
oe (7) 
T(ktat+lF(k+jtat+6 +1) ; 
= 2 Fite ATG lad Drees py be 
Thereby, we deduce the following: 
(a,8);p, _ Tk+a+1) 
Prk (L) — T(a +1) k! ’ 
‘papi. UP Tet Bt het et B+): 
BYP ce O) DT(k—-r+)F@+B+l) (8) 
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- TikK+atlj(kK+a+6+4+1), 
De = 10 
Pai el L(k—-r+1)0(r+atl) ’ mM) 
rp(@,8) = T(r 4p Oe B a5 1) atr,B+r) 
D = . 11 
Prk (x) LT(kta+B+1) Pin (x) (11) 
Taking wor? ) (x) = (L—2)°2*, we list the following norm and inner product 


related to the weighted space: L* vo.) (0, £] as 
L£ 


L 
(1,0) fo = f we) te w(x) dz, |lvll yea = (0,0)? non (12) 
0 


A complete L? (a,s) (0, £J-orthogonal system is consisted of a set of shifted 
Ww 


L£ 
Jacobi polynomials, where 


ENP ei) lool 
Pe als we 8) — “@) hy é = he k . (13) 
We used x4 a a: ) anda oe eae ) 0 <7 <N, as the nodes and Christoffel numbers 


of the standard Jocabes Gauss eae ate in the interval [—1, 1]. 
The corresponding nodes and corresponding Christoffel numbers of the 


shifted Jacobi-Gauss interpolation in the interval [0, £] can be given by 


Lyte 
_ ~@ +1), 


le 
TLNj ~ 9 


,B) 
Nj 
a L a a,B : 
Deny = (GPM wy S, OS GSN. 


For any positive integer N, @ € Soy41[0, £] and by means of Jacobi-Gauss 
quadrature property, we obtain 


fe £—2)*2° ba) yar = (5) ° a] 


(—2)*(1+2)%¢ (Fe ) dx 


o+B8+1 N 
L a L (a, 
z =) ore (Fey i 1)) (14) 
j=0 
~ (a, 8) (a, 8) 
= Vet? (xf in 
j=0 


Use Caputo’s fractional derivative given in (3) and (4) for the shifted 


Jacobi polynomials POG t) to obtain its fractional derivative as 


pDEP&-(t)=0, 7=0,1,...,fu]—-1, w>0, (15) 
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0,84) Ly (aa a, 814) ~n— 7 
Po Dapp PE => gloPMen—w 5 = Tu], fu] +1,..., (16) 


n=[H] 


where c #) ig given by 
Tg+64+)DF(n+j+at84+1)t 


T(int+ B+ UVPGtat+Bt+ ly —-—n)'P(n-pt lr 


CS ( 1)-” 


3 One-dimensional T-FOTE 


3.1 One-dimensional T-FOTE with homogeneous 
conditions 


Here, we are occupied with utilizing the SJC method to solve one-dimensional 
T-FOE [29]: 


O'u(a, t) 


Ov-1u(z, t) 
r YL 
OtH 


Ote-1 


+ yu(x,t) = 73 (17) 


where 1 < yp < 2 and the term H(2,t) denotes the field variable. We also 
assume the homogeneous conditions: 


(0,4) = u(£,t) = u(e,0) = MEP), <0, (18) 


OMu(a, t) 
OtH 


derivative. In addition, 71, y2 and y3 are given constant coefficients. The 
point of our method is to get solution can be extended, using combination of 
basis functions of shifted Jacobi polynomials, in the form 


where £ and 7 are given and represents the Caputo fractional 


N-2 M-—2 
u(x,t) = D> Sc, po'2*? (ypler) (4) = (a, tC, (19) 
i=0 j=0 


where, ¢,;, 7 =0,1,...,N—2, 7 =0,1,...,M — 2 are the unknown coeffi- 
cients, 


C = [co,0, €o,1, -- - > C0, M—23C1,0; C1,1)-- C1, M—23 


7 
CN—2,0,CN—-2,1)--> ,CN-2,M-2| , 
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N and M are any arbitrary positive integers, and 


Gr FY (op) = PLP (a) 4 eg PLLAY (x) + eg PLOY (x), (20) 
EPO) = PLP O + PLEO + aPRRPO. 2D 


From the boundary conditions ¢ eas (0) = perPu(c) = 0 and the two 

relations (8), we have the accompanying framework 
_ti+)) -, @ti+ G+ 1+ 2) - 
“ G4+1) °° G@4+1)G42) 


i, (22) 


tit), 41+ YG+142) _ 
* G41) °° G@4+1)G4+2) 


Thus €; and €; can be remarkably resolved to give 


- PNGU=1(e+141 
G+14+1G+14+D(Qi+14 
(i (i+ 2)(+141 
G+14+1G@4+14+)D(Qi+1 


(24) 


y= 


Also, one can without much of a stretch check that 


247 + 1)(27 +2+2+4+3) 

29 +2+2+ 4)’ 

j+2)(27 +2424 2) 
g+242)(27 424244) 


(25) 


Fy 
GHC 
or Ge 21 


From (15) and (16), we have that 


Q2,P2, a2, 82 2,82, 2;P2, 
org = Dire Oa Pre Otero 
+ oe), 


and 


a2, 82,— _ a2,B2 a2, 82,6 (a2,82,u— 
ie b De) = Di tpl i ie t) =p Bayt Vie t) + pyPie2 Pe: 7 Dee) 


(a2,82,4—1 
+ oyPoatorD ey, 


(27) 
Also y(a,t) is the 1 x (N — 1)(M — 1) matrix introduced as follows: 


(a, t) =|roo(a,t), 70,1, t),< +5 70,at-2(@, 8); r1,0(&, t), 71,1(2, 2), 


: .,71,M—2(2,t); TN—2,0(2,t), rw—2,1(2, t), te ,TN-2,M—2(2,t)], 


where 
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rig (a, t) = GSEP) (aw) pOr")(t), 4 =0,1,.. 


Substituting (19) into (17) yields 


gu N=2M=2 


a Ye cadet”? wre”) 


i=0 j=0 


_1 N-2M-2 
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aN=2 4=0.1 =e 


tngger( DD aud” er) 
i=0 j=0 
N-2M~2 5) 
O1;P1 (a2, 82) 
a er, Tied aC 
i=0 j=0 
N-2 M-2 
=o ‘em De: Ci, j ey bie pier Pa) (t )) + H(e, 2), 
1=0 j=0 
N-—2 M-—2 
O1;P1 or a2 , 2: 
Cig ot (@) sa a) 
7=0 g=0 
N-2M-2 aun 
a a2, 2 
+m cag 6P (@) ples 
i=0 j=0 
N-2M~2 (29) 
hp Cie? OWeeP'O) 
i=0 j=0 


(@1,81) 
Cig Lyi 


form: 


(x) br") (t)) + H(a, t). 


Therefore, adopting (26) and (27), enable one to write (29) in the following 
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N-2M-2 
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S- Ss Ci.5 eg :P1) (x) per Bo, Ht t) 


i=0 g=0 


N-2M-2 


1, aa a 
+n >> S- cis¢et 1,81) ( (x) pe 2825p De) 


i=0 g=0 
N-2M-2 


(30) 
$72 > YS es 9h3 wy (t) 


1=0 7=0 


‘>b>) 


nF) a ypienP®) (t)) + H(a, t). 
1=0 j=0 
Assume that 


fig (2, t) = GP a ypenPary ) +44 ae on, B1) Guerre 


+ nder? (eRe) - wo mele ) 


noe); 
at that point, (30) can be modified as 


N-2M-—2 


yi ci fig (@, t) = 
7=0 7=0 


U(x, 1). (31) 


Collocating (31) in N—1 and M—1 roots of the shifted Jacobi polynomials 
Pe ,(z), the Gauss-shifted Jacobi nodes, we obtain 


a1,81 a2,B2 a1,81 a2, Bo 
Cy fig (2 vi ) fee ¥ ee — H(x ( ) th ) 
7=0 jg=0 


nyt IT mg ); 


forn=0,1,...,N—2, m=0,1, 


d 
which can be written in the following matrix form 


FTC=R, 
where 


= [Roo, Rio,---, Rn—2,0;Ro1, Ria,---,Rn—2,13 


Ro,m—2; Ri,m—2;---; Rn—2,M—-2]" 
Rig = H(oeuty), he"), ¢=0,1,...,N—2, f= 0,1,...,M—2, 
and 


F = (fijnm); t,n=0,1,...,N—2, j,m=0,1,...,M—2 


bs 
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in which the elements of the matrix F are determined as follows: 


fijnm = fis: an — DS ncng iN ey j,M= Oy Necng AE = 2, 


In our implementation, this system has been solved using the Mathemat- 
ica function FindRoot with zero initial approximation. In this manner, the 
approximate solution of (17) is given by u(z,t) = y(a,t)C. 


3.2 One-dimensional T-FOTE with nonhomogeneous 
conditions 


In the accompanying, we simplify alter the right-hand side to deal with the 
nonhomogeneous initial-boundary conditions. Give us a chance to consider, 
for example, one-dimensional time-fractional order telegraph equation (17) 
with the nonhomogeneous initial-boundary conditions: 


u(0, t) =qo(t), u(L, t) = q(t), O0<t<T, 


Ou(a, t) (33) 


u(x, 0) =po(x), at lt=0 = pi(2), O<a2< £, 


where qo, 91, Po; P1 are known functions and the function u is unknown. 
Presently, assume the accompanying transformation 


V(a,t) = u(a,t) + ao(t) + ay (t)x + a(x — L£L)(bo(x) + bi (x)t), (34) 
where 


ag(t) = —go(t), ar(t) = go(t) — a(t) 


L ? 
so(a) = L=2)M0(0) + (0) = Lol) 
. La(x —L) : 
_ 2a lo — Pi (x) # (2ag(4) lio Sa) r=) 
es as a(a —L) ; 


The mapping (34) changes the nonhomogeneous conditions (33) into the 
following homogeneous conditions: 


V(0,t) = V(L,t) = V(2,0) 


Subsequently it suffices to solve the following time-fractional order tele- 
graph equation: 
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OKV(a,t) | OP! V (a, t) — PVE) 
ate get + BV (at) =~ — a + G(@,8): 


O<e<£,0<¢<T, 


(36) 


with the homogeneous conditions (35), where V(x, t) is obtained by (34) and 


a) arte 4) (Oraolt) , . O#*anlE 
G(x, t) =H(x,t) + ( aie aes 


O#ay(t) , Otay (t) 
id nr 


+ a(x — £)(b1 (x) re) T3-p)’ 


T(2—») 
+ a0(lo(e) + bx(0)t) — yo( APO? 4 PHO) 


Mota) + M1(2) 273(bo (x) + b1(x)t))). 


) + Y2a9(t)) 


T(2 
t'—# 4 aby (2) ( ) 2-p 


— 2y3(2x — £)( 


4 Two-dimensional T-FOTE 


4.1 Two-dimensional T-FOTE with homogeneous 
conditions 


In this section, we test the following two-dimensional T-FOTE: 


du(a,y,t) _ O¥tu(z, y, t) 


Ott r V1 Otel y2u(zx, y, t) 
Pu(a,y,t)  Pu(zy,t), , (37) 
= y3( Ox2 I Oy? ) t Hie, Yy, i), 
(2e<f. <9 <0 <'<7, 
with the homogeneous conditions 
u(0,y,t) = u(Li,y,t) = ula, 0,t) = u(x, Lo, t) = u(x, y, 0) 
_ Ou(z, y, t) (38) 


y=, 
at Ji=0 


The point of our method is to get solution can be extended, using combination 
of basis functions of shifted Jacobi polynomials, in the form 
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where, ¢;,;,4, 7=0,1,...,N — 
the unknown coefficients, 


C= [co,0,0, CO,0,1)-+- 
€1,0,0,C1,0,1,--- 
CN —2,0,0;CN—2,0,1)+-+5 


T 
-,CN—2,M—2,K—2]" 5 


2,j7=0,1,... 


»€0,0,K —2,C0,1,0)C0,1,1;-++ 
»€1,0,K —2,€1,1,0,€1,1,1;--+ 


CN—2,0,K —2;CN—2,1,0)CN—2,1,1)--- 
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(a) pS? (ype?) (2) 


»C0,1,K—2)+++;C0,M—2,K—23 


»C1,1,K—2)-+-,C1,M—2,K—2; 


»CN—2,1,K—2) 


N, M, and K are any arbitrary positive integers, and 


a1,h1 a1,R1 a1,P1 a, 81) 
OEP? (a) = PLP? (a) + GPL (0) + eePEHY (@), (40) 
a2, a2, a2, a2,8 
gor wy ) = Poa (y ) = Pea (y Nag ej, Poy ), (41) 
Ure (t) = Pre) + eePreer (1) + oP pees (t), (42) 


where €, €;, €;, €j, Pk and g, are defined in (24) and (25). 


Also x(x, y,t) is the 1 x (N — 1)(M — 1)(K — 1) matrix introduced as 
follows: 

x(a, Y, t) =(ro,0,0(z, Y) t), r0,0,1(2, Y, t), Cae! »70,0,K—2(2, Y, t), ro,1,0(2, Y, t), 
ro,1,1(2, Y, t), adie ,T0,M—2,K—2(2, Y, t); T1,0,0 (x, y, t), 
POET eye) ole ait) 

15,4 Ub) ye234 Ti Roo We) Pao eG) 
rn—2,0,0(£, Y, t), Tn—2,0,1(2, yy; t), PRS TN—2,0,K—2(2, Yy, t), 
rn—2,1,0(2, Y, t), rn—2,1,1(2, Y, t), ee) TN—2,1,K—2(2, yy, t), 994 
TN-2,M—2,K-2(2,y, t)], 
where 
; 3,83 : 
neu =e  @ee yO) tee Oh t=O Ga —2, 
j=0,1,....M@—2, k=0,1,...,K —-2. 


Substituting (39) into (37) yields 
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N-2 M-2 K-2 
a1, Pr) ae al (a3, 3) 
Sa Do Ye esate (aelers eyed? p) 
i=0 j=0 k=0 
juni N=2M-2K~2 


Sat eevee ; 63 Ss 5 ape (x) (@2,82) 


a3,83 

se )pho2Pa) (yap PP) (zy) 
i=0 jg=0 k=0 

N-2 M-2 K-2 


a1,P1 2 ,P2 03,83 
$72 > SY cp no'?? a) eer (yee (0) 


i=0 j=0 k=0 


(43) 
g2 Na2M=2K-2 


=e a o> S> yo ae et (x) (a2 ,B2) 


x) pie) (yyy lor) (4)) 
i=0 j=0 k=0 
a? N-2M-2 K-2 


1,P1 2,P2 a3,83 
1 ait (YY sco @eery wre) 


i=0 j=0 k=0 
+ H(z, y,t), 


2 
5 
w 
> 
w 


Cj Pre) ie ae (x) po2sP) ae ai yes 83) (t) 


ll 

i) 
a. 

ll 
° 


a 


| 
3 
z 
w 
= 
w 
f 


a1, a2, ov 03,83) 
cj ePene? aoe (y Y) aaa pies PN) 


l 
Co 
& 
wn S 
n> 


Co 


z 
g 


nN 


| 
T 
No 
Ke) 
& 
> 
as 
ee 


cee ae (Were? O 


© 
Il 
oO 
GS 
Il 
oO 
> 
Il 
o 


(44) 


® 
re) 
: 
x 
w 
x 

wo 


° 


3 
= 
bd 
vu 
a 


enon” Ooo wwe oO) 


Q 

Nw 
2t 

| 


1;P1 a2,P2 a3,83 
+55 ( Cin @)PES WYER O) 
Y" 7=0 j=0 b=0 


+ H(z, y,t). 


Therefore, adopting (40), ( 


form 


41), and (42), enables one to write (44) in the 
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2 
I) 
Sd 
w 
x 
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c-2 


Ye ck Oee me wre PO) 
k= 


L2,J 
0 


Il 
i) 
&. 
ll 
° 


2 
f 
=< 
a 
“A 


a1, a2, B: a3,B3,u—1 
+71 on er” (x ep e222) (yap! ed } )(t) 


Co 


yk 


i) 


a1,P1 2 ,P2 03,83 
+ Cig ROE PY (a) plE2P (ype? (4) 


° 
co 
ll 
fos) 


(45) 


® 
i) 
i 
x 
. 
x 

E 


2 
S 
ak 
M 


Cin Oene (een) We eE” ®) 


°. 


@ 
Nw 
2b 


1,P1 2,P2 a3,83 
+7855 Cig RPE EY (a) per (ywree™ (t)) 


Assume that 


fae(®, 9st) =O0 2? (@) pene? Were O) 


+ nde eer? Were YO 


$y \022P) (ar) yo'022P2) (yy gp loP) (t) 
02 


8 Oa a (Oe meee? OvrR ?O) 


L2,J 


a a2, a3,63 
Nine SSE wreSeP wwe”), 
at that point, (45) can be modified as: 


N-2M-2 K-2 


PS », » GH ahem th=H(a, Yy, t). 


(46) 
i=0 j=0 k=0 


Collocating (46) in N —1, M—1 and K — 1 roots of the shifted Jacobi 


polynomials eo (x), the Gauss-shifted Jacobi nodes, we obtain 


7d (a1, Bi) (a2, B2) (a3, 83) 
Cij.kbigk(@e, 7 ni Yeo, mJ try k ) 
i=0 j=0 k=0 


(47) 
=e ea ), for n=0,1,...,N—2 
m=0,1,...,M—-—2, l=0,1,...,K -2, 


which can be written in the matrix form F7C = R, where 
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R = [Ro,0,0, Ro,o,1,---, Ro,0,K-2) Ro,1,0; Ro,1,1,-+-,-Ro1,K-2) +++; 


Ro,m—2,K—2; 21,0,0, R1,0,1,---, Ri0,K—-2, 11,0, Ri1j1,---, 

Ri ijK-2,---,Ri,m—2,K-23 Rn-2,0,0, Rn—-2,0,1,---, RN—2,0,K-2; 
T 

Ry-2,1,0, Rn-2,1,1,---,RN-2,1,K-2,---,RN-2,M—2,K-2]° , 


Q1, a2, B: a3,B3 . 
Rete ae a), Ca ha 
j=0,1,...,M—2, k=0,1,...,K —2, 
and 


F =(fijknmi), 1,n=0,1,...,N —2, j,m =0,1,...,M —2, 
a One ee 


in which the elements of the matrix F are determined as follows: 


: a2, 3,63 
fijknml =e. 2 oo "), 


in=0,1,...,N—2, 
j,m=0,1,...,M—2, k,l =0,1,...,K —2. 


In our implementation, this system has been solved using the Mathemat- 
ica function FindRoot with zero initial approximation. In this manner, the 
approximate solution of (37) is given by u(z,y,t) = x(z,y,t)C. 


4.2 Two-dimensional T-FOTE with nonhomogeneous 
conditions 


In the accompanying, we simplify alter the right-hand side to deal with the 
nonhomogeneous conditions. Give us a chance to treat with two-dimensional 
(2D) T-FOTE (37) with the nonhomogencous initial-boundary conditions: 


u(0,y,t) =go(y,t), wLi,y,E=al(y,t), OS y<L2,0<t<T, 
u(x, 0,t) =qo(z,t), u(x, £2,t) =q3(z,t), O< ae <L1,0<t<T, 


Ou(a, y, t) 
at 


(48) 
lt=-0 =Q5(@,y), OS @<SL1,0<y< Lo, 


u(x, y, 0) =q4(z,y), 


where qo(y; t), 71 (y, t), q2 (x, t), 93 (x, t), qa(z, y) and q5 (a, y) are given func- 
tions. 


Presently, assume the accompanying transformation 

V(z,y,t) =u(x,y,t) + ao(y,t) + rar (y, t) 
+ &(2 — £1) (do(x, t) + ybi (2, t)) (49) 
“Ir x(x _ Li)yly _ L)(co(x, y) ale te, (2, y)), 
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where 
qo(y,t) — a(y,t 
ao(y,t) =—doly,t), arly, t) = - ), 
Li - 0,t)-—L£ t 0.t 
bo(a, t) = (Li — £)qo(0,t) — Ligo(z, t) + 7m (0, ) 
a(x —Ly)Ly 
bs (w, t) = Lilge@t) — a(@,#)) + (Lr — #)(go(La,t) ~ 0(0, 4) + 2 (Lat) ~ 1,4) 
~ LiLox(x — Li) ; 
co(%, ¥) 


L£1(q2(, 0) + go(y, 0) — aa(z, y)) + x(qr(y, 0) — a1(0, 9) — go(y, 0)) — 0 (0, 0) (L1 — 2) 
ryLly(x — L1)(y — £2) 

, YEi(93(@, 0) — ga(#,0)) + y(L£i — ©)(go(0, 0) — go(L2, 0)) + ey(qi (0, 0) — a1 (£2, 0)) 

wvyliLo(a — Li)(y — £2) : 


ei (2, y) 


jee ae tesa) = 849 (0,#) Saale. t) } a( 24 (y 8) _ Sao (yt) 840 (0.t) = ae 


ayLly(« — £4)(y — £2) -_ 


pS — 29206) 4 y(cy — 2)( 20000) _ Pan(E 2) 4 py OM O.4) _ | 
t=0 


wylyLo(@—L£1)(y — £2) 


The mapping (49) changes the nonhomogeneous conditions (48) into the 
following homogeneous conditions: 


V(0,y,t) =V(Li,y, t) = V(a,0,t) = V(a, Lo, t) = V(2, y, 0) 
_ VEE) a (50) 
= at t=0 — V; 
Subsequently it suffices to solve the following 2D T-FOTE: 
OVvVia,y,t) , OF OV (a, 9, #) 
Ott rY1 otu-} y y2V (a, y, t) 
_ (FV wt) _ PV (2, y,t) 
ae Ox? Oy? 
O0<2<L,,0<y<Ll2,0<t<T, 


)+G(x,y,t), OY 


subject to the homogeneous initial-boundary conditions (50), where V(z, y, t) 
is given by (49), and 
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OK a9 (y, t) 
Ott 


dH tao(y, t) 
inat + 1240(Y,#) 


— 273 (bo(a, t) + ybi(z,t)) — 73 Fooly, t) qo a) 


Oy? ote 
Ol ai(y, t) 0? a1(y, t) 
Otu-1 


OM bo (a, t) Ob; (a, t) 
+a(e —L£i)( Ait ner 
O#~1bo (a, t) O#—1b, (a, t) 
Otu-1 Ote-1 ) 
+ Y2(bo(x, t) + ybi (2, t)) 
07 bo (x, t 07by (x,t)  d%co(a, 
rr of, 1 / Vy olf wy ie 
Ox Ox Oy Oy 
dco(x, y) dei (x,y) 
+ 2 2t 


ale — Ly )(y = Le) (eta, y) 


G(x, Y, t) =H(a, Y, t) + 


+1 + y2a1(y, t) — 3 


+ y1( 


r(2)je!-# 
P(2 — p) 
Fe + valealtsy) + ter(2yy) 
ee) 

— 245 (22 f)( Poles | ytd) 

— 273 (y — £2) (co(z, y) + ter(z, y)) 


— Qy3(2x — Li) (y £)( Aor) 


+ 1¢1(2,y) 


Oci (2, y) 
t Aa ). 
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5 Convergence and error analysis 


In this section, following the analysis given in [11], we have the following 
theorems. It is worthy to note here that the proofs of the next theorems are 
similar to the proofs given in [11]. 


Lemma 1. The Caputo fractional derivative of the shifted Jacobi poly- 
nomials satisfies the following estimate: 


CD" PEP (@)| seats, 


where C is a positive generic constant and g = max(a, , —3). 
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(26 +21 + Do 
Lita, +1)¢+6i+1) 


, then we have the following two connection formulas: 


Theorem 1. If €; = 


(27 +A2+ De 
= 727 + Bo + l)a 


Gilt) = a(L — a) PLETED (e), 


w; (t) =a PP rn (x). 


and 


Proof. The two formulas are easily obtained by application of the moment 
formula given in Doha [7]. 


Theorem 2. The following two orthogonality relations are valid 


L 
[ bi(x)b;(2) 2-1 (L — 2) de 


— £63 (24+A1 41)? (24+A14+2)TGt+a4+)DPG+64+1) 
_ (@+o, +1) (4+ 6,+1) TG 4 14 2) 


? 


[ " in (t)n (t) 192-2 ( — 1)? at 


LA2-? (2m + Az + 1)? (2m + Az + 2)P(m+ ag +: 1)P(m+ Bo +1) 
(m+ Bz + 1)2m!T(m + AQ + 2) 


= dmn 


Proof. The proof is a direct consequence from the orthogonality relation (13). 


Theorem 3. (Convergence) If u(x,t) the exact solution of (17) is separa- 
ble in the sense that u(x,t) = 2(L—ax)tf(x)g(t) and f,g are C® functions 
with | f’ (x) |< a, |g’ (t) |< b, where a,b are positive constants, then the 
expansion coefficients c,;; in (19) satisfy the following estimate: 


| ci; |= O(i 2972) for all i,j > 3. 
Moreover the series in (19) converges absolutely as N, M — oo. 


Theorem 4. If u,un,y are the exact and approximate solutions of (17), 
respectively, and under the assumptions of Theorem 3, then we have the 
following truncation error estimate 


llu — un,mll2 = O(N? M2). 


Theorem 5. If en, = u—uy,m is the truncation error of the solution of 
(17) and under the assumptions of Theorem 3, then we have the following 
global error estimate: 
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| DY en. +1 °DE* en.m +72 eNn,m — 73 D2 ew alo 
<p, M¥t9-2 N~2 + pa |yi|Met9-2 N-? 
+p3 ly2| M~2 N~2 + pa|q3| M72 NOt? 


Theorem 6. [Convergence] If u(x,t) the exact solution of (37) is separable 
in the sense that u(x, y,t) = t(Lo — x)y(Le — y)t f(x) g(y) A(t) and f,g,h 
are C® functions with | f(x) |< a, |g (x) |< b, | h(t) |< c, where a,b,c 
are positive constants, then the expansion coefficients ¢;;, in (39) satisfy the 
following estimate: 


| cija |= O(@- 2-2-2) for all i,j,k > 3. 
Moreover the series in (39) converges absolutely as N, M, K > ov. 


Theorem 7. If u,un,w,x« are the exact and approximate solutions of (37), 
respectively, and under the assumptions of Theorem 6, then we have the 
following truncation error estimate: 


Iu — unm, llo = O(N7~2M-2 K-2). 


Theorem 8. If en wjK = U-—UN,m,xK is the truncation error of the so- 
lution of (37) and under the assumptions of Theorem 6, then we have the 
following global error estimate: 


| (DY enue +1 °DE* en.m.x +72 eNn.M,K — 73 (D2 + D*) en,.M,K 
< py KO+9-3 M-2 N-? + fy|y|Kett-? M-2 No? 

+s \y2|M~? N~? K~? 

+hsl (p12 M~2 Net? K-43 + fy N-3 etd K-?) 
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6 Numerical results 


In this section, a few examples are outlined to demonstrate the pertinence 
and proficiency of the novel technique in homogeneous and nonhomogeneous 
conditions. The calculations are executed by utilizing Mathematica of Ver- 
sion 8, and all counts are completed in a PC of CPU Intel(R) Core(TM) 
i3-2350M 2 Duo CPU 2.30 GHz, 6.00 GB of RAM. 

The distinction between the measured value of approximate solution and 
its actual value (absolute error), are given by 


E(z,t) = |u(zx, t) _ u(z,t)|, (52) 
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E(az,y,t) = |u(z,y,t) — a(a, y, £)|, 


(53) 


where u(x,t) and u(x,t) are the exact solution and the numerical solution 
at the point (x,t), respectively. Also u(z,y,t) and u(z,y,t) are the exact 


solution and the numerical solution at the point (z, y,t), respectively. 


Moreover, the maximum absolute errors (MAEs) is given by 


MAEs = Max{E(z,t) : V(z,t) € [0, £] x [0, 7]} = L™, 


MAEs = Max{FE(a, y,t) : V(x, y,t) € [0, £1] x [0, £2] x (0, 7]} = L. 


Also we can denote to the root mean square error (RMSE) by 


N-2 
> (ulaler®), 2) — talon) 2)2 
RMSE = \| ° 


N 
or 
N-2M-—2 
B a2 By ~( (O18 1B 
a p> wee es t) ~~ ‘ae? ’ Tees t))? 
RMSE = \| —— 


NM 


(54) 


(55) 


(56) 


(57) 


Example 1. Give us initial a chance to consider the accompanying T-FOTE 
(17) in the domain 0 < x < £, with taking after initial and boundary condi- 


tions, 


u(0,t) =u(£,t)=0, O<a<, 
_ Ou(z, t) 
Ot 


u(a, 0) lio = 0, t>0. 


The exact solution is given by 
u(x,t) = t?(£ — x)sin?(z), 
and 


2t?-# (2 — £L) sin? (x) (u — yt — 3) 
P(4— p) 
+ 2t?73((a — L) cos(2x) + sin(2z)) 


H(x,t) = 


t to (L — 2) sin*(x) 


RMSE is shown in Table 1 for different values of a1, 81, a2, G2, V1, Y2 
and y3 and H(x,t). Figure 1 shows the exact solution and numerical solutions 
of Example 1, where N = 14, M=4, 9 =72 =73=1, L=T=5, pw 
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Table 1: Comparison of RMSE for Example 1, where £ = JT = 1 and N = 12, M =4. 


V1) 72573 Our method (2 = 1.25) Our method (j: = 1.75) 
ay ~ 1 : y ay a = ay = 1 ~ ft ay =e = 
a2 = Be 3 a2 = Pa= 5 a2 = By 2 a2 = Be 3 


71 = 0,72 =73 = 1] 5.8845.10-" | 2.6905.10-1° | 6.0047.10-"% | 8.1272.10-" 
y2 =0,71 =73 = 1 | 1.7662.10-' | 4.6050.10-'° | 6.3574.10-' | 1.2146.10~ 1° 
73 =0,% =72 =1)| 1.8565.10-" | 2.4882.10-! | 9.0080.10-! | 7.2685.10-! 
a =1,4%2=73=0)] 1.7330.10-" | 2.2806.10-" | 1.0194.10-" | 1.1746.10-" 
92 = 1.4% =73 =0 |] 2.2531.10-" | 2.3049.10-" | 6.4857.10- | 8.7682.10-” 
=ly=%=0/ 3.3257.10-" | 2.6122.10-% | 8.29299.10-" | 2.6548. 10-% 


Figure 1: The exact u(z,t) and numerical U(z,t) solutions for Example 1, where N = 


144,M=4, 1 =72 =73 =1, C= T =5, and wp = 1.95. 


1.95, and a, By ag = Bo 0. In Figures 2 and 3, we depicted the 
following: 


Fig. 2 The curves-graph of exact u(x,t) and numerical u(x,t) solutions for 
Example 1, where N = 14, M=4, y =72=73=1, L=T7=5, pw 
1.95 and a, = 8; = ag = G2 = 0 at five different values of t, 


Fig. 3 The curves-graph of exact u(x,t) and numerical u(x,t) solutions for 
Example 1, where N = 14, M =4, y, =72 =73 =1, L=T=5, ps 
1.95 and a, = 8, = a2 = Pz = 0 at five different values of zx. 


It should be noted here that the CPU time 7 used to implement our 
algorithm was computed using the TimeUsed command in Mathematica and 
was found to be 15 < 7 < 22, 17 < 7 < 28 in seconds for N = 12, M = 4, 
N = 14, M = 4, respectively. 

Example 2. In this example, we considered the T-FOTE (17) with the 
coefficients 7, = y2 = 1, y3 = am and the following initial and boundary 
conditions: 
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/ 
wx.) ¥ 0.5.1) 
iix,1) f 0.5.1) 
3 Wx,2) a ul.5.t) 
2 ~ f¥x.2) 3 f ~ is 
3 | ) / is, 
. x3) ft w.5.) 
~ ii3) of - 302.53) 
7 Pa 
ix,4) if vG.5,t) 
iix,4) 3.5.1) 
wix,5) “ ae ws) 
~ ffx.) ttt ~ (4.5.1) 
Lo | 1 2 3 4 5 


Figure 2: The curves-graph of exact u(x,t) and numerical u(x,t) solutions for Example 
1, where N= 14, M=4, 71 = 72 = 73 = 1, L=T =5, and uw = 1.95 at five different 
values of t. 


Figure 3: The curves-graph of exact u(x,t) and numerical %(2, t) solutions for Example 
1 where N = 14, M=4, 71 72 = 73 1, £=T7 =5 and p = 1.95 at five different 
values of x. 


u(0,t)=0, ulLoi=fsm (Ll), 0<¢<7, 


lr-0 = 0, O<a2<. 


where the source term H(z, t) is chosen to be consistent with the exact 
solution u(a, t) = #3 sin?(z). 

The results of Example 2 are reported in Table 2 and Figures 4, 5, and 6. 
Table 2 compares RMSE with [26] at different values of uw, N, M, ai, 81, a2, 
and (2. The space-time graphs of the exact and approximate solutions at 
N = 20, M = 3, £=T = 10, and pw = 1.05 are plotted in Figure 4. We 
plotted in Figures 5 and 6. 


Fig. 5 The curves-graph of exact u(x,t) and numerical u(x,t) solutions for 
Example 2, where N = 20, M = 3, £=T7 = 10, and wp = 1.05 at five 
different values of t, 


Fig. 6 The curves-graph of exact u(x,t) and numerical u(x,t) solutions for 
Example 2, where N = 20, M = 3, £=7 = 10 and pw = 1.05 at five 
different values of x. 


It should be noted here that the CPU time 7 used to implement our 
algorithm was computed using the TimeUsed command in Mathematica and 
was found to be 33 < 7 < 41 for N = 20, M =3. 

Example 3. Let us test the T-FOTE (17) with the coefficients 7, = y2 = 
y3 = 1, and H(z,t) according to the exact solution u(x,t) = xcos(x? + t?) 
and the following initial and boundary conditions: 
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Table 2: RMSE for Example 2 at L= 7 =1. 


Lb Sweilam et al. Our method (NV = 12, M =3) 
(m = 20, n=3) [a1 =~ =F, | 1 =f, =0 on = Bi = 3, 
[26] a2 = Pp =—5 | A2=f2=0 | ag = fo =5 
1.75 1.0474.10~° 4.2800.10-" | 7.4066.10-™ | 7.8044.10~-™ 
1.95 1.0080.10~° 1.8994.10-!! | 2.7463.10-'! | 9.4081.10~" 
m Sweilam et al. Our method (VN = 12, M =4) 
(m = 20, n = 4) a a, = 6; =0 on = Bi = 3, 
[26] a2 = B2a=— 3 | 2 =f2=0 | op =f. =35 
1.75 1.0429.10~° 8.2456.10- 1 | 7.9494.10-™ | 2.2020.10~ 1 
1.95 9.8685.10~° 3.5018.10-1! |] 1.3740.10-!° | 3.8586.10~ 18 


Figure 4: The exact u(x,t) and numerical u(2,¢) solutions for Example 2, where N = 
20, M=3, L=T = 10, and p = 1.05. 


umd 1 
u and fi 


Figure 5: The curves-graph of exact u(x,t) and numerical «(z, t) solutions for Example 
2, where N = 20, M =3, £L=T = 10, and p = 1.05 at five different values of t. 


Figure 6: The curves-graph of exact u(x,t) and numerical u(z, t) solutions for Example 
2, where N = 20, M=3, L=T = 10, and up = 1.05 at five different values of x. 
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Table 3: TRMSE for Example 3 at £= 7 =1. 
js | Hosseini et al. Our method (N = 12 = M = 3) 

(N=n=50) | a, = 8, =—, | a =f, =0 an = Bi = 5, 

[18] a2 = Bg =—5 | 2 = fo =0 | a2 =fPo= 5 

1.25 | 2.6551.1073 5.7636.10-8 | 5.9213.10-® | 5.8662.10—§ 
1.5 4.1743.10-% 3.4483.10~§ 5.8363.10-5 | 1.0047.10~7 
1.75 | 5.4164.10-3 4.2018.10-8 9.7424.10-8 | 2.2745.10-7 
1.95 | 3.8800.10-% 1.4632.10- | 4.6185.10-® | 4.3415.10-8 


Figure 7: The exact u(x,t) and numerical w(z,¢) solutions for Example 3, where N = 
20, M = 20, £=T7 =3, and wp = 1.05. 


u(Oi)=0, u(Lji)=Len(l? +e), Off <7, 
u(x,0) = xcos(z”), Pate =0, O0<a<L 


In Table 3, we give the RMSE with various choices of 4, a1, 61, a2, 
Bo, and £ = JT = 1. The outcomes are contrasted with the outcome of 
radial basis function [18]. It is clear from this table that the solution gotten 
by our technique is great in examination with radial basis function [18]. The 
numerical and exact solutions are compared in Figure 7, where N = M = 20, 
a, = By = ao = 62 = 1, w = 1.05, and £ = T = 3 Moreover, in Figures 8 
and 9, we sketched the following: 


Fig. 8 The curves-graph of exact u(x,t) and numerical u(x,t) solutions for 
Example 3, where N = M = 20, L=T =3, a) = By = a2 = fo = 1, 
and uz = 1.05 at five different values of t, 


9 The curves-graph of exact u(x,t) and numerical u(z,t) solutions for 
Example 3, where N = M = 20, L=T =3, ay = 0} =a2 = fo = 1, 
and ys = 1.05 at five different values of x. 


Fig. 
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Figure 8: The curves-graph of exact u(x,t) and numerical u(x,t) solutions for Example 
3, where N = 20, M = 20, £=7 = 3, and up = 1.05 at five different values of t. 


Figure 9: The curves-graph of exact u(x,t) and numerical %(z, t) solutions for Example 


3, where N = 20, M = 20, £=7 =3, and p = 1.05 at five different values of x. 
Table 4: RMSE for Example 4 at pp = 1.5. 

1=1=-3, 1 =1= 0, 1=1= > 

(N, M, k) 2=2= —-39; 2=2= 0, 2=2= 3 

ee =.= 0 eg 1 

3 3 3 3 3 3 =3= 5 
(2,22) | 4.8245x10-> | 1.2683x 10-7 | 1.9055 x 10-2 
(4,4,4) | 2.3555 x 10-5 | 2.1954 x 10-5 | 2.0665 x 10-5 
(6,6,6) | 1.5245 x 10-7 | 1.6220 x 10-7 | 3.9732 x 10-7 
(8,8,8) | 2.9602 x 10-!° | 3.4822 x 10-2¢ | 3.9732 x 107-22 
(10,10,10) | 3.7346 x 10-13 | 5.1405 x 10-18 | 7.0877 x 10-8 


Example 4. We consider the 2D T-FOTE (37) with the coefficients 7, = 
Jo = 73 = 1, and H(z, y,t) according to the analytical solution u(z,y,t) = 
t? sin(x + y) and the following initial and boundary conditions: 


u(0,y,t) = t? sin(y), u(1,y,t) =t?sin(y+1), O<y<1, 0<t<1, 
u(x,0,t) = ¢?sin(x), u(a,1,t) =tsin(e +1), O<2<1,0<t<1, 


du(a, y, t) 
Ot 


u(x, y,0) = 0, L-9 = 0, O< a <1, O<y<1. 


Tables 4 and 5 show the RMSE and the MAEs, respectively, with different 
values of a, 3, N, M, and kK. 
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Table 5: MAEs for Example 4 at 4 = 1.5. 

1=1=-3, 1 =1= 0, 1=1= > 
(N, M, Kk) 2=2= 5, 2 =2= 0, 2=2= 9 
a er) 3 =3= 0 3 =3> 5 

(2:90) 4.0454 x 10-8 | 3.7132 x 10-3 | 3.4097 x 1073 

(4,4,4) 4.1167 x 10-° | 3.7308 x 10-° | 3.4409 x 107° 

(6,6,6) 2.1882 x 10-7 | 2.2566 x 107” | 2.2929 x 1077 

(8,8,8) 4.5656 x 107-!° | 5.5665 x 1071 | 6.4892 x 10710 

(10,10,10) | 5.6932 x 10713 | 9.0227 x 10718 | 1.3389 x 10718 


7 Conclusions 


The SJC method is effectively connected for finding the agreement of T- 
FOTEs. We have accomplished a decent understanding between the approx- 
imate solution acquired by SJC and the exact one. The aftereffects of the 
examples demonstrate that the shifted Jacobi collocation method is depend- 
able and effective technique for solving time fractional telegraph equation and 
also other equations. The key feature of the proposed method is to obtain 
highly accurate semi-analytic solutions via few number of retained modes. 
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